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SECTION  1 


INTRODUCTION 

A  high  altitude  nuclear  explosion  exposes  the  upper  atmo¬ 
sphere  to  fast  photoelectrons  produced  by  absorption  of  the 
x-rays  from  the  bomb,  and  to  beta  rays  (megavolt  electrons) 
from  the  radioactive  debris.  These  electrons,  and  the  second¬ 
ary  electrons  they  produce,  directly  dissociate  some  of  the 
oxygen  in  the  air,  and  also  create  ions  and  other  reactive 
species  that  lead  to  further  dissociation.  The  atomic  oxygen 
thus  created  is  of  practical  importance  because  it  detaches 
electrons  from  negative  ions,  which  increases  the  absorption 
of  radio  waves.  In  addition,  it  combines  with  molecular 
oxygen  to  form  ozone  (Og),  which  is  a  strong  infrared  emitter, 
particularly  in  the  9  to  11  fim  spectral  region. 

Current  predictions  of  the  amounts  of  atomic  oxygen  and  ozone 
produced  per  unit  energy  deposited  in  the  upper  atmosphere  are 
based  on  rather  complex  calculations  of  the  relevant  electron- 
impact  and  chemical  processes,  since  no  direct  measurements 
are  available.  However,  several  laboratory  measurements  of 
the  ozone  production  in  irradiated  air  and  oxygen  at  pressures 
typical  of  the  lower  atmosphere  have  been  carried  out  and 
published.  It  is  useful  to  utilize  these  experimental  results 
to  check  or  improve  the  complex  air  irradiation  calculations , 
especially  since  the  calculations  involve  several  uncertain 
cross  sections  and  rate  constants.  A  logical  procedure  is  to 
treat  the  pure  oxygen  case  first,  without  the  complication  of 
nitrogen  or  water  vapor  reactions.  This  is  done  in  the  pres¬ 
ent  analysis;  even  without  the  other  air  components  it  yields 
some  significant  results  pertinent  to  atmospheric  situations. 

A  later  report  will  treat  the  normal  air  case. 
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A  few  previous  researchers  have  made  limited  comparisons  of 
oxygen  radiolysis  measurements  with  chemical  calculations  . 

The  most  recent  and  most  complete  of  these  is  that  published 
by  Willis  et  al .  (Ref.  1)  in  1970.  Although  these  investi¬ 
gators  were  able  to  obtain  reasonable  agreement  of  calcula¬ 
tions  with  measurements,  they  assumed  rather  than  calculated 
the  initial  number  of  oxygen  atoms  produced  by  the  irradia¬ 
tion.  Moreover,  although  they  used  the  best  chemical  rate 
coefficients  available  at  that  time,  some  of  their  values  are 
now  known  to  be  inaccurate,  and  they  neglected  excited  state 
reactions.  Accordingly,  a  new  and  more  complete  analysis  is 
needed  in  order  to  confirm  or  correct  the  conclusions  of 
Ref.  l. 
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SECTION  2 


INITIAL  SPECIES  PRODUCTION  IN  THE  RADIOLYSIS  OF  OXYGEN 

2 . 1  COMPUTATIONAL  PROCEDURE . 

The  first  step  in  making  ozone  calculations  for  comparison 
with  the  measurements  is  to  determine  the  species  produced 
initially  by  the  electron  bombardment.  This  requires  a  com¬ 
plex  calculation  of  the  excitation,  dissociation  and  ioniza¬ 
tion  processes  that  occur  as  the  primary  and  secondary  elec¬ 
trons  slow  down  to  thermal  energies.  Several  groups  of 
researchers  have  developed  computer  codes  for  making  such 
calculations  for  any  gas  of  interest,  provided  that  the  rele¬ 
vant  electron-impact  cross  sections  are  input.  Two  groups 
have  published  a  few  results  for  pure  oxygen  (Refs.  2,3). 
However,  neither  group  has  given  a  set  of  results  complete 
enough  for  the  present  purposes,  nor  did  either  employ  the 
best  available  input  cross  sections. 

Accordingly,  new  calculations  of  the  initial  species  produc¬ 
tion  in  oxygen  have  been  made  at  RDA.  The  results,  are  sum¬ 
marized  in  Table  1 .  The  RDA  calculations  employed  a  "discrete 
energy  bin"  code,  DPOSIT,  which  is  based  on  a  code  originally 
developed  by  Dalgarno,  Lejeune  and  Victor  (see  Ref.  4),  modi¬ 
fied  by  Stephens  and  Klein  to  handle  higher  energy  primaries 
(Ref.  5),  and  further  modified  for  greater  input  and  output 
flexibility  by  the  present  writers  and  coworkers  at  RDA.  It 
has  been  checked  against  a  comparable  code  developed  by  Green 
and  colleagues  (Ref.  2)  by  applying  it  to  irradiated  air,  and 
found  to  give  results  agreeing  within  2  percent  in  every 
detail,  when  the  same  input  cross  sections  are  used. 
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Table  1.  Initial  production  of  atoms,  ions  and 
excited  states  in  bombarded  oxygen*. 


Species 

G  Value  ( particles/ 100  eV) 

O2U1  V 

8 . 4  * 

°2(blEg) 

1.7* 

0(3P) 

4.0 

0(  *D) 

3.0 

0+2 

2 . 25 

0+ 

1 .0 

e 

3 . 25 

*  The  results  are  essentially  identical  for  bombardment  by 
1  keV  and  10  keV  primary  electrons. 


t  These  values,  which  were  used  in  the  present  chemistry 
calculations,  were  computed  with  0.15*  CO2  *n  'the 
02.  In  pure  oxygen  the  values  would  be  slightly  larger, 
8.9  and  1.8,  respectively,  but  this  makes  a  negligible 
difference  in  the  ozone  chemistry. 
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2.2  IMPORTANCE  OF  ELECTRON- IMPACT  CROSS  SECTIONS. 


The  principal  source  of  error  in  electron  degradation  calcu¬ 
lations  is  the  inaccuracy  in  the  basic  electron-impact  cross 
sections.  Using  limited  measurements  and  theory,  Porter, 
Jackman  and  Green  (Ref.  6)  derived  algebraic  expressions 
for  the  electronic  excitation  cross  sections  of  oxygen, 
and  Jackson,  Garvey  and  Green  (Ref.  7)  did  similarly  for 
the  ionization  cross  sections.  However,  more  recent  measure¬ 
ments,  as  well  as  some  older  measurements  which  they  over¬ 
looked,  show  that  their  expressions  are  frequently  signifi¬ 
cantly  in  error,  sometimes  by  more  than  a  factor  of  two.  For 
the  present  analysis,  a  careful  search  for  and  evaluation  of 
the  best  available  cross  sections  was  made,  and  the  results 
were  input  to  the  DPOSIT  code  as  numerical  tables,  avoiding 
the  approximation  of  algebraic  fits . 


A  complete  DPOSIT  calculation  requires  inputs  and  provides 
outputs  for  a  large  number  of  excited  and  ionized  states  . 

For  present  purposes,  the  most  important  outputs  are  the 
total  amounts  of  oxygen  atoms  and  0+  and  0^  ions  produced. 

In  addition,  the  low  metastable  states  of  02  and  0  are 
involved  in  ozone-destroying  reactions,  and  so  should  be 
included.  Moreover,  several  other  excited  states  also  need 
to  be  included  in  the  DPOSIT  calculation,  because  excitation 
removes  energy  from  the  colliding  electrons,  leaving  less 
energy  for  dissociation  or  ionization  by  subsequent  colli¬ 
sions.  In  addition,  dissociation  of  molecules  by  electron 
impact  actually  occurs  by  excitation  of  unstable  electronic 
states  which  spontaneously  dissociate;  inclusion  of  such 
states  is,  of  course,  essential  to  the  present  analysis. 
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Since  so  much  depends  on  input  cross  sections,  the  sources  of 
the  data  used  in  the  present  calculations  will  be  summarized 
below  in  some  detail,  starting  with  the  lowest  electronic 
states  and  working  upward. 

2.3  CROSS  SECTIONS  FOR  NEUTRAL  EXCITATION  AND  DISSOCIATION. 

For  the  low  02  metastable  a1^  and  b^E  states,  smooth  curves 
were  drawn  through  the  data  of  Linder  and  Schmidt  (Ref.  8)  , 
Trajmar  et  al .  (Refs.  9  and  10),  and  Wakiya  (Ref.  11),  which 
are  generally  consistent  in  their  overlapping  energy  regions . 
For  the  combined  A,  C  and  c  states  near  5  eV ,  which  are  too 
close  together  to  be  resolved  in  electron-impact  measurements , 
the  data  of  Refs.  10  and  11  were  used.  Since  no  measurements 
below  15  eV  incident  energy  are  available,  the  cross  section 
was  assumed  to  fall  rapidly  below  this  energy,  in  analogy  to 
the  rapid  decrease  of  the  cross  section  for  the  neighboring 
B  state,  discussed  below. 

A  curve  for  the  important  B3£  state,  which  contributes  over 
half  of  the  total  oxygen  atom  production,  was  drawn  through 
the  data  of  Trajmar  et  al.  (Ref.  10)  and  Wakija  (Ref.  12)  . 

A  similar  curve  for  the  low  Rydberg  states  between  9.7  and 
12.1  eV  was  based  on  the  data  of  Wakiya  (Ref.  12).  For  both 
states  the  published  data  extend  only  down  to  incident  ener¬ 
gies  of  20  eV .  Below  this  energy  both  cross  sections  were 
assumed  to  decrease  rapidly,  in  accordance  with  a  preliminary 
unpublished  value  by  Shyn  and  Sharp  (Ref.  13)  for  the  B  state 
at  15  eV . 

As  mentioned  earlier,  electron-impact  dissociation  of 
molecules  generally  takes  place  by  excitation  to  unstable 
electronic  states  which  spontaneously  dissociate.  In  oxygen 
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this  is  believed  to  occur  for  essentially  every  excitation 
to  the  A  +  C  +  c  and  B  states,  and  to  every  Rydberg  state 
below  12.1  eV  (the  ionization  energy  of  02).  For  the  A  +  C 
+  c  and  B  states  this  conclusion  follows  from  the  Franck- 
Condon  principle  applied  to  the  potential  curves  (Ref.  14) 
and  also  from  the  observed  electron  energy  loss  (Ref.  10-12) 
both  of  which  indicate  production  of  these  states  above 
their  dissociation  limits,  so  that  the  two  oxygen  atoms 
immediately  fly  apart.  For  the  Rydberg  states,  the  observed 
broadening  of  the  ultraviolet  absorption  lines  and  the  failure 
to  detect  corresponding  emission  lines  indicate  that  these 
states  "predissociate ,  11  i.e.,  spontaneously  dissociate  by 
interaction  with  neighboring  dissociating  states.  The  cross 
sections  for  all  of  these  dissociating  states  could  be  com¬ 
bined  into  a  single  "dissociation  cross  section."  However, 
because  the  energy  lost  by  the  incident  electron  varies  with 
the  particular  state  excited  (the  excess  over  the  dissociation 
energy  going  into  kinetic  or  excitation  energy  of  the  atomic 
fragments),  it  is  desirable  in  DPOSIT  calculations  to  treat 
each  dissociating  state  separately. 

In  addition  to  the  states  discussed  above,  Porter  et  al  . 

(Ref.  6)  also  considered  a  number  of  higher  Rydberg  states. 
Since  none  of  these  states  are  observed  in  optical  emission, 
they  must  either  predissociate  or  autoionize.  Porter  et  al  . 
arbitrarily  assumed  a  50-50  split  between  these  two  processes 
for  most  of  the  states,  while  photoabsorption  and  photoioniza¬ 
tion  data  (Ref.  15)  show  that  autoionization  is  generally 
favored.  In  addition,  the  estimates  of  Porter  et  al .  for  the 
cross  sections  of  these  states  were  based  on  their  estimates 
of  the  cross  sections  for  the  excited  states  of  0£ 
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(the  Rydberg  "cores"),  which  they  generally  overestimated,  as 
discussed  below.  Consequently,  their  Rydberg  cross  sections 
are  probably  too  high. 

It  would  be  difficult  to  re-estimate  cross  sections  for 
all  of  the  18  higher  Rydberg  states  treated  by  Porter  et 
al . ,  and  the  results  would  still  have  considerable  uncer¬ 
tainty.  Fortunately,  these  states  contribute  less  than  5 
percent  to  the  total  dissociation  cross  section  and  less  than 
15  percent  to  the  total  ionization  cross  section.  Moreover, 
the  ionization  error  can  be  largely  removed  by  constraining 
the  total  of  all  the  ionization  cross  sections  to  the  measured 
total  ionization  value.  Accordingly,  for  the  present  DPOSIT 
calculations  only  the  two  lowest  (and  most  strongly  excited) 
Rydberg  states  above  12.1  eV ,  i.e.,  the  (a4n)3J?  and  (A2n)3J? 
states  at  about  12.7  and  13.3  eV,  were  included.  The  cross 
sections  used  were  those  Porter  et  al .  (Ref.  6).  Although, 
as  discussed  above,  these  cross  sections  are  probably  too  large, 
this  is  offset  by  the  neglect  of  higher  Rydberg  states.  In 
accordance  with  photoabsorption  and  photoionization  data 
(Ref.  15),  these  states  were  assumed  to  dissociate  spontan¬ 
eously  30  percent  of  the  time  and  autoionize  70  percent  of 
the  time . 

2.4  CROSS  SECTIONS  FOR  IONIZATION. 

The  cross  section  for  ionizing  O2  to  its  ground  ion  state, 

X  2n  ,  was  derived  by  subtracting  the  cross  sections  for 
Rydberg  autoionization,  just  discussed,  and  for  excited  ion 
states,  discussed  below,  from  the  total  ionization  cross  sec¬ 
tion.  The  total  cross  section  was  taken  from  the  measurements 
of  Rapp  and  Englander-Golden  (Ref.  16)  up  to  200  eV  and  of 
Schram  et  al .  (Refs.  17,  18)  above  500  eV ,  with  a  smooth 
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transition  in  the  intermediate  energy  region  where  the  data 
of  the  latter  lie  about  10  percent  lower  than  those  of  the 
former . 

Turner  et  al .  (Ref.  19)  measured  the  fraction  of  metastables 
in  an  0+  beam  produced  by  bombarding  02  by  electrons  of  ener¬ 
gies  up  to  200  eV.  This  fraction  rises  from  zero  near  16  eV 
to  a  value  of  0.33  above  70  eV .  These  results  can  be  con¬ 
verted  to  a  cross  section  for  producing  metastables  by  multi¬ 
plying  by  the  total  oxygen  ionization  cross  section,  deter¬ 
mined  above,  and  by  the  fraction  of  the  total  ions  ( 0 +  +  o+ ) 
that  are  0  4 ,  measured  by  Rapp  et  al .  (Ref.  20).  The  only 
known  or  predicted  metastable  state  of  0+  is  the  a4n  state  at 
16.11  eV  above  the  ground  state  of  02  (Ref.  14).  This  0^ 
state  can  be  produced  either  directly  or  by  the  rapid  radia¬ 
tive  decay  of  the  b4!  and  c4E  states.  However,  the  c-state 
apparently  contributes  little  to  the  metastable  production, 
since  it  primarily  dissociates  (Ref.  6).  Accordingly,  the 
a-state  cross  section  has  been  obtained  by  subtracting  from 
the  metastable  cross  section,  the  cross  section  for  producing 
the  b-state,  derived  below.  It  might  be  noted  that  this  new 
a-state  cross  section  is  about  two-thirds  of  that  employed 
by  previous  investigators  (Refs.  2,  3,  5-7). 

The  cross  sections  for  exciting  a  few  emission  bands  of  the 
0+  A2n  state  have  been  measured  by  Korol'  et  al .  (Ref.  21), 
while  a  single  point  for  one  band  at  100  eV  has  been  given  by 
McConkey  and  Woolsey  (Ref.  22).  However,  the  ratios  of  seve¬ 
ral  band  intensities  given  by  Korol1  et  al .  disagree  seriously 
with  the  ratios  obtained  in  radiative  studies  (Ref.  23) . 
Moreover,  the  cross  section  of  Korol'  et  al .  falls  off  below 
100  eV  with  unrealistic  steepness.  Fortunately,  both  this 
data  and  that  of  McConkey  and  Woolsey  show  that  the  A-state 
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cross  section  is  less  than  one  percent  of  the  total  ionization 
cross  section,  and  so  can  be  neglected  for  present  purposes. 
(The  correct  cross  section  is  clearly  over  an  order  of  magni¬ 
tude  smaller  than  the  semiempirical  cross  section  derived  in 
Ref .  7)  . 

The  cross  section  for  exciting  the  0+  b4E“  state  has  been 

2  g 

measured  with  modern  methods  by  five  different  groups  of 
investigators  (Refs.  22,  24-27),  omitting  the  discordant 
measurements  of  Korol'  et  al .  (Ref.  21)  for  reasons  discussed 
earlier.  All  obtained  curves  of  essentially  the  same  shape, 
with  a  maximum  around  100  eV,  but  their  absolute  magnitudes 
are  significantly  different,  presumably  due  to  differences  in 
normalization  (often  a  problem  in  this  field).  After  correc¬ 
tion  for  incomplete  coverage  of  all  the  b-state  emission 
bands,  the  cross  sections  of  three  of  the  investigating  groups 
(Ref.  22,  24,  25)  agree  within  10  percent,  while  that  of  Borst 
and  Zipf  (Ref.  26)  is  about  40  percent  lower,  and  that  of 
Srivastava  (Ref.  27)  is  about  50  percent  lower. 

As  discussed  earlier,  the  sum  of  the  a-  and  b-state  cross 
sections  is  fixed  by  the  metastable  measurements  of  Turner 
et  al .  (Ref.  19).  If  the  b-state  cross  section  of  Refs. 

22,  24  or  25  is  assumed  correct,  the  two  cross  sections  are 
found  to  be  roughly  equal,  while  if  that  of  Borst  and  Zipf 
(Ref.  26)  or  Srivastava  (Ref.  27)  is  used,  the  b-state  cross 
section  is  found  to  be  about  1/2  or  1/3  that  of  the  a-state . 

Guidance  in  choosing  among  these  alternatives  can  be  obtained 
from  optical  measurements.  Theory  and  measurements  (Ref.  28) 
show  that  above  a  few  hundred  eV  the  electron-impact  cross 
section  for  ionization  is  approximately  proportional  to  the 
photoionization  cross  section  (or  differential  optical 
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f-number)  integrated  over  photon  energies  with  a  weighting 
factor  that  emphasizes  the  lower  energies.  Samson  et  al . 

(Ref.  29)  have  measured  the  partial  photoionization  cross 
sections  for  producing  the  b-state  and  the  combined  a  +  A 
states  (unresolved),  between  threshold  and  41  eV .  Over  most 
of  this  region  the  a  +  A  cross  section  is  about  twice  that 
of  the  b-state;  moreover,  it  is  quite  large  in  the  16-18  eV 
region,  which  is  below  the  threshold  of  the  b-state.  Since, 
as  discussed  above,  optical  emission  data  show  that  the  con¬ 
tribution  of  the  A  state  to  the  a  +  A  combination  is  small  , 
one  concludes  that  the  electron-impact  cross  section  for  pro¬ 
ducing  the  a-state  is  at  least  twice  that  for  producing  the 
b-state.  This  rules  out  the  b-state  cross  sections  of  Refs. 

22,  24  and  25.  Of  the  two  remaining  sets  of  data,  that  of 
Borst  and  Zipf  (Ref.  26)  has  been  chosen,  because  they  are 
more  experienced  workers  and  their  data  are  smoother  . 

(These  values  are  only  about  half  as  large  as  the  semiem- 
pirical  values  of  Ref.  7. ) 

Electron  impact  on  0 also  produces  several  excited  states 
of  0+  above  the  b-state.  Apparently  all  of  these  states 
spontaneously  dissociate  into  0+  +  0  (Ref  6),  except  for  a 
small  fraction  of  the  c-state,  which  can  be  neglected  for 
present  purposes.  Hence,  they  are  included  in  the  dissocia¬ 
tive  ionization  cross  section  measured  by  Rapp  et  al .  (Ref.  20) 
Much  less  is  known  about  the  partition  of  this  cross  section 
among  the  different  excited  states,  but  fortunately  the  DPOSIT 
calculations  are  less  sensitive  to  the  partition  than  to  the 
total.  Jackman  et  al .  (Ref.  7)  used  the  limited  experimental 
evidence  to  divide  these  states  into  three  groups,  with  thres¬ 
holds  at  20,  23  and  37  eV .  Their  groups  and  cross  sections 
have  been  used  in  the  present  analysis,  except  that  some 
adjustments  were  made  to  the  20  eV  cross  section  to  make 


11 


the  total  of  the  three  agree  with  the  measurements  of  Rapp 
et  al .  ( Ref .  20 ) . 

2 . 5  RESULTS  OF  ELECTRON  IMPACT  CALCULATIONS . 

The  results  of  the  DPOSIT  calculation  are  summarized  in  table 
1,  with  the  following  modifications:  All  predissociating 
excited  states  (discussed  above)  are  replaced  by  their  disso¬ 
ciation  products,  0  +  0  or  0+  +  o.  The  ground  state  (3P)  and 
metastable  ( )  atoms  are  tallied  separately,  which  requires 
information  about  the  electronic  excitation  of  the  dissocia¬ 
tion  products.  The  lowest  dissociating  states,  A  +  C  +  c,  and 
the  dissociatively  ionizing  states  below  20.7  eV,  have  insuf¬ 
ficient  energy  to  produce  excited  0  atoms.  The  well-studied 

^2  B  state  is  known  to  dissociate  into  0(3P)  +  0(*D).  The 
products  of  the  higher  excited  states  are  less  well  estab¬ 
lished,  but  they  contribute  only  about  20  percent  of  the  total 
0  atoms.  Consequently,  for  these  states  the  product  assign¬ 
ments  of  Porter  et  al .  (Ref.  6)  have  been  used  without  further 
study.  The  small  amount  of  0(1S)  product  they  suggest  is 
assumed  here  to  be  quenched  immediately  to  0(3P)  . 

The  state  of  excitation  of  the  ionized  products,  0+  and  0+ , 
has  been  ignored  for  present  purposes,  since  on  collision  with 
°2,  excited  02  rapidly  quenched  and  excited  0+  rapidly 
transfers  its  charge  (Ref.  30).  Such  reactions  might  possibly 

yield  metastable  0^  or  o.  but  since  there  is  no  evidence  for 
this  possibility,  it  has  been  neglected.  The  DPOSIT  calcula¬ 
tion  gives  G(0+)  =  2.42  and  G(0  +  )  =  1.04,  where  the  G  values 
are  defined  as  the  number  of  particles  produced  per  100  eV 
of  deposited  energy.  The  sum  of  these  values,  3.46,  slightly 
exceeds  the  experimental  value,  G(ions)  =  3.25,  derived  from 
several  careful  measurements  in  oxygen  showing  that  30.8  eV  is 
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required  to  produce  each  ion  pair  (Ref.  31).  Although  the 
difference  is  only  6  percent,  the  experimental  value  appears 
accurate  to  about  1  percent.  Accordingly,  for  table  1  the 
calculated  values  of  G(0+)  and  G(0+)  have  been  decreased  by 
6  percent  to  give  agreement  with  the  ionization  measurements . 

The  DPOSIT  calculation  was  carried  out  for  1  keV  primary 
electrons,  while  the  oxygen  radiolysis  experiments  generally 
used  bombardment  energies  of  hundreds  of  keV.  However,  theo¬ 
retical  calculations  (Ref.  6)  and  ionization  measurements 
(Refs.  31,  32)  indicate  that  the  G  values  are  constant  within 
about  5  percent  for  primary  energies  above  1  keV . 
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SECTION  3 


CHEMICAL  REACTIONS  IN  OXYGEN  RADIOLYSIS 


3.1  SOURCES  OF  BASIC  DATA. 

The  chemical  reactions  believed  to  be  involved  in  oxygen 
radiolysis  are  listed  in  table  2,  together  with  their  rate 
coefficients.  Much  of  the  chemical  kinetic  data  is  taken 
from  standard  compilations  of  laboratory  measurements  (Refs. 
33-38),  and  some  from  more  specialized  papers  (Ref.  39-49)  . 
Several  rate  coefficients,  however,  have  not  been  measured, 
so  have  had  to  be  estimated  by  analogy  with  similar  reactions 
or  by  use  of  indirect  evidence,  as  indicated  in  the  footnotes 
to  the  table . 

3 . 2  ION-MOLECULE  PROCESSES . 

In  oxygen  radiolysis,  the  initial  bombardment  produces  0+ ,  0+ , 
and  thermal  electrons.  The  electrons  very  rapidly  attach  to  0^ 


The  primary  0+  ions  also  rapidly  charge  exchange  with  0^: 


(2) 


0+  +  0 


+  0  . 


The  0“  and  0+  ions  then  evolve  by  charge  transfer,  ion-neutral 
clustering,  and  other  processes.  The  basic  clustering  reac¬ 
tions  involve  0^; 


(3) 


+  o2  -  o4  +  o2 


(4) 
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Table  2.  Reactions  and  Rate  Coefficients 
for  the  Pure  Oxygen  System. 


Reaction 

Rate  Coefficient^3' 

Reference  or  Footnote 

e  +  0*  +  0  +  0(1D) 

2.1  x  10'7(300)0,63 

T 

33 

0++  02  +  o2  +  o 

1.9  X  10"11 

34 

02  +  02  +  02  +  O4  +  02 

4.0  X  10'30  (30°)  2,93 

T 

39 

O4  +  O2  *  O2  +  2^2 

1.3  X  10"4  (300)2. 93^-6000/T 

b 

e  +  0*  +  0  +  0  +  02 

2.0  X  10"6  (300) 

33 

oj  +  0  +  02  +  03 

3.0  X  10"10 

34 

0*  +  O3  +  02  *  0^  +  02 

1.0  X  10"28  (300)3 

c 

°5  +  °2  *  °2  +  °3  +  °2 

3.3  X  10'4  (300)3  - 7850/ T 

T  e 

d 

e  +  0^  >  0  +  202 

2.0  X  10'6  (30°) 

T 

e 

°4  +  °3  *  °5  +  °2 

1.0  X  10'10 

f 

°5  +  °2  *  °4  +  °3 

1.0  X  10'11  e_1850/T 

f 

e  +  02  +  02  >  02  +  02 

1.4  X  10"29  (300)  -600/T 

T  e 

33 

02  +  0  +  e  +  03 

1.5  X  10'10 

33 

02  +  0  >  0  +  0^ 

1.5  X  10'10 

33 

°2  +  °2  (XAg)  +  e  +  202 

2.0  X  10_1° 

34 
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Table  2.  Reactions  and  Rate  Coefficients 
for  the  Pure  Oxygen  System  (Continued). 


Reaction 

Rate  Coefficient 

Reference  or  Footnote 

0"  +  0  e  +  0^ 

2.0  X  10"10 

33 

°’+  °2  "  e  +  °3 

3.0  X  10"10 

34 

°2  +  03  +  03  +  02 

O 

X 

1 — * 

O 

1 

1 — * 

O 

33 

03  +  0  y  0^  +  ^2 

2.5  X  lO"10 

34 

O3  +  O3  y  O2  +  20^ 

3.0  X  io-1V1000/T 

g 

0^  +  0^  +  ^  ^4  +  ^2 

3.5  X  10-31  (300) 

33 

0^+0  y  O3  +  0^ 

4.0  X  10-10 

34 

O4  +  O3  y  O3  +  203 

O 

r-H 

1 

O 

rH 

X 

O 

CO 

33 

0+  +  0"  >  0  +  02 

2.0  X  10"6 

h 

0++  0^  >  0  +  202 

2.0  X  10-6 

h 

02  +  0">  0  +  02 

<sD 

1 

O 

t— H 

X 

0 

C\J 

h,i 

O2  +  O2  +  0  +  0  +  0^ 

2.0  X  10"6 

h 

°2  +  °3  *  °2  +  °3 

2.0  X  10-6 

h 

+  — 

O2  +  O4  **  0  +  0  +20^ 

ro 

0 

X 

1 — » 

0 

1 

or 

h 

oj  +  0'  +  0  +  202 

2.0  X  10-6 

h 
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Table  2.  Reactions  and  Rate  Coefficients 
for  the  Pure  Oxygen  System  (Concluded). 


Reaction 

Rate  Coefficient 

Reference  or  Footnote 

Oj  +  O2  ■*  0  +  0  +  202 

2.0  X  10"6 

h 

O4  +  O3  *  03  +  202 

2.0  X  10-6 

h 

O4  +  O4  >  0  +  0  +  3O2 

2.0  X  10-6 

h 

0^  +  +  2O2 

2.0  X  10"6 

h 

0*  +  O3  +  2O3  +  o2 

2.0  X  10'6 

h 

0^  +  O4  ^  +  3O2 

2.0  X  10"6 

h 

0  +  02  +  02  *  O3  +  02 

6.9  X  10"34  (300)1-25 

35 

0  +  O3  +  02  +  02 

8.0  X  10-12  e"2060/T 

36 

o  +  o  +  o2  >  o2  +  02 

4.8  X  10'33  (300)2 

j 

0  (XD)  +  02  *  0  +  02  O-Z*  ) 

2.6  X  10-11  e67/T 

36 

0  (*0)  +  02  >  0  +  o2 

6.0  X  10"12  e67/T 

37 

0  (XD)  +  O3  ■*  02  +  02 

1.2  X  10"10 

36 

0  (XD)  +  O3  +  0  +  0  +  02 

1.2  X  10-10 

36 

02  (^g)  +  02  *  02  +  02 

2.2  X  IQ'18  (-4)0-8 

38 

0  .l  +  l 

u2  (  \)  +  02  +  02  (  Ag)  +  02 

3.8  X  10"17 

47,48 

02  (lA  )  +  03  +  0  +  202 

1.2  X  lO-1^-240071 

38 

U2  (  Zg)  +  03  >  0  +  202 

2.4  X  10"11 

49 
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Footnotes  to  Table  2 


a.  Rate  constant  units  for  biraolecular  reactions  are  cra3/sec, 
and  for  termolecular  reactions,  cm6/sec. 

b.  The  reverse  rate  coefficient  is  calculated  using  the 
forward  rate  coefficient  and  the  equilibrium  constant  data 
of  Payzant  et  al .  (Ref.  40)  and  Conway  and  Janik  (Ref.  41), 
Keq  =  3.0  x  10-26  e600°/T  cm3. 

c.  The  value  for  He  as  the  third  body  is  about  1  x  10-28 
cm6/sec  at  200  K  (Ref.  34).  The  same  value  is  used  for 
02  as  the  third  body  at  300  K.  Further,  a  temperature 
dependence  of  (300/T)3  is  assumed,  as  in  the  02  +  02  +  02 
reaction  (Ref.  39). 

d.  Calculated  from  the  reverse  rate  coefficient  and  the 

equilibrium  constant  derived  from  the  equilibria  for  02  + 

°2  +  °2  *  +  °2  an<*  +  °3  *  +  ®2  ( see  footnotes  b 

and  f )  . 

e.  Assumed  to  have  the  same  value  as  for  e  +  0^[ . 

f.  The  forward  reaction  is  estimated  from  other  fast  ion- 
neutral  switching  reactions  (Ref.  34);  the  reverse  reaction 
is  calculated  from  the  equilibrium  constant  of  Dotan  et  al . 
(Ref.  42),  Keq  =  10  e185(5/T. 

g.  This  rate  coefficient  was  selected  to  have  a  room  temper¬ 
ature  value  of  about  10-12  cm3/sec,  as  deduced  in  oxygen 
radiolysis  experiments  (Ref.  43),  and  an  activation 
energy  of  about  2  kcal/mole,  as  deduced  from  temperature- 
dependent  studies  (Ref.  44).  The  present  analysis  sug¬ 
gests  that  this  value  may  be  an  upper  limit,  and  a  more 
likely  value  is  about  half  as  large  (see  section  4.4). 

h.  Estimated  high  pressure  value  from  References  45  and  46. 

i.  It  is  possible  that  the  products  are  0+0+0,  but  this 
reaction  appears  to  be  unimportant  both  in  radiolysis 
experiments  and  in  the  ambient  or  disturbed  atmospheres. 

j.  The  rate  coefficient  listed  corresponds  to  N2  as  the  third 
body  (Refs.  33  and  38).  Insufficient  data  are  available  at 
temperatures  below  500  K  with  02  as  the  collision  partner. 


19 


The  importance  of  the  clusters  is  principally  in  their  greater 
thermochemical  stability,  which  might  prevent  them  from  break¬ 
ing  up  to  yield  oxygen  atoms  when  they  recombine  with  oppo¬ 
sitely-charged  ions  (see  next  section).  Ion  molecule  inter¬ 
change  reactions  can  also  lead  to  ions  of  greater  stability: 

°4  +  °3  "*  °5  +  °2  •  (5) 

Such  reactions  have  not  been  considered  in  past  radiolysis 
studies . 

In  the  case  of  negative  ions,  molecules  with  large  electron 
affinities  can  scavenge  negative  charge  by  charge  transfer. 

Two  such  reactions  that  are  important  in  oxygen  radiolysis  are 


°2  +  °3  "*  °3  +  °2  ' 

(6) 

and 

°4  +  °3  °3  +  °2  +  °2- 

(7) 

Previous  radiolysis  analyses  have  omitted  the  latter  reaction. 

Certain  reactions  can  detach  electrons  or  reduce  complex  ions 
to  simpler  ones,  for  example: 


°2  +  0  ->  e  +  03  , 

(8) 

CM 

O 

+ 

1  CM 

O 

t 

O 

+ 

1  00 

0 

(9) 

o4  +  0  -*  o3  +  o2  . 

(10) 

However,  at  the  high  pressures  typical  of  radiolysis  experi¬ 
ments,  such  recycling  reactions  are  overwhelmed  by  clustering 
reactions . 
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At  pressures  near  one  atmosphere,  larger  cluster  ions  such 
as  ot  and  o7  are  known  to  form,  but  few  of  their  reaction 

D  D 

rate  coefficients  are  known.  Equilibrium  data  (Ref.  41) 
indicate  that  their  concentrations  are  small  except  below 
200°K ,  so  they  have  not  been  included  in  the  present  analysis. 

3.3  ION-ION  RECOMBINATION. 

In  radiolysis  experiments,  the  charge  produced  by  irradiation 
is  neutralized  primarily  by  recombination  of  positive  with 
negative  ions.  In  the  present  case,  the  positive  ions  are 
predominantly  02  and  0* .  and  the  recombination  products  can 
be  02  or,  if  enough  energy  is  available  for  dissociation,  0+0 
The  distinction  is  important  because  the  0  atoms  generated  by 
ion-ion  recombination  can  eventually  form  ozone  and  increase 
the  ozone  yield.  If  all  recombination  occurs  via  stable  ions 
without  dissociation,  no  net  ozone  production  from  ions  results 

In  the  most  recently  published  discussion  of  the  ozone  yield 
under  various  conditions  of  oxygen  radiolysis  (Ref.  1),  it  was 
suggested  that  reactions  of  the  following  sort  lead  to  0-atom 
generation : 


(°2  or  °t)  +  (°2  °r  “4)  ' 


0  +  0  +  n  02 , 


(ID 


while  the  following  reactions  do  not: 


^°2  or  °4 1  +  03  ■*  03  +  n  °2  •  (12) 

These  suggestions  have  been  adopted  in  the  present  calcula¬ 
tions.  We  also  assume  that  the  recombination  of  ot  with 

o 

negative  ions  leads  to  02  +  03  without  0-atom  production. 
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Measurements  of  ion-ion  recombination  at  atmospheric  pressure 
suggest  a  typical  recombination  coefficient  of  about  2  x  10-6 
cm^/sec  (Refs.  45,  46).  A  specific  study  has  been  made  by 
Wood  et  al .  (Ref.  50)  of  ion-ion  recombination  rates  in  oxygen 
subject  to  low  level  continuous  gamma  irradiation.  At  pres¬ 
sures  of  100-300  torr,  an  average  coefficient  of  about  1.1  x 

_C  a 

10  cm°/sec  was  measured.  At  a  pressure  around  30  torr  the 
coefficient  decreased  to  about  0.9  x  10-6  cm3/sec.  Since  most 
radiolysis  experiments  have  been  conducted  at  pressures  near 
one  atmosphere,  in  the  present  calculations  a  single  ion-ion 
recombination  coefficient  of  2  x  10-fi  cm3/sec  has  been  adopted 
for  all  pairs  of  ions  at  all  pressures  and  temperatures. 

3.4  OZONE  CHEMISTRY. 


The  principal  neutral  reactions  controlling  the  ozone  concen¬ 
tration  are  the  production  mechanism: 


0  +  °2  +  02 


°3  +  32  ' 


and  the  destruction  mechanism: 


(13) 


0  +  °3  -  02  +  °2 •  (14) 

Such  reactions  can  occur  for  both  ground  and  excited-state 
oxygen  atoms,  with  different  rate  coefficients.  In  the  radi- 
olysis  experiments  of  interest,  excited  oxygen  atoms  are  so 
rapidly  quenched  by  oxygen  molecules  that  only  the  ground- 
state  reactions  need  to  be  considered. 
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Ozone  can  also  be  destroyed  by  ion-molecule  reactions,  such  as 


°2  +  °3  -  Oj  *  02,  (6| 

°i  *  °3  -  °2  *  202-  <15) 

This  reaction  sequence  is  "catalytic*1  in  that  the  second  reac¬ 
tion  regenerates  the  0“  lost  in  the  first  reaction.  A  similar 
pair  of  destruction  reactions  could  involve  0  +  and  0*  in  place 
of  0”  and  0”,  respectively. 

The  excited  molecular  oxygen  states  02(1a)  and  09 ( )  can 
dissociate  ozone: 


02  (  1A)  +  °3  ■*  0  +  2  02  ,  (16> 

o2( h  )  +  03  -  0  +  2  o2.  (IV) 

However,  at  the  pressures  typical  of  radiolysis  experiments, 
almost  all  of  the  oxygen  atoms  thus  produced  quickly  react 
to  re-form  ozone.  Hence,  the  role  of  these  excited  states 
is  minor . 

An  upper  limit  to  the  amount  of  ozone  that  can  be  produced  by 
radiolysis  may  be  derived  by  neglecting  all  ozone  destruction 
processes.  With  this  neglect,  every  oxygen  atom  that  is 
created  yields  one  ozone  molecule.  From  table  1,  for  each 
100  eV  of  energy  deposited  in  0^  7.0  oxygen  atoms  are  pro¬ 
duced,  and  also  1.0  0+  ions,  which  rapidly  charge  exchange 

with  02  -j-0  yield  another  oxygen  atom.  Hence,  from  these 
processes  alone,  G(O^)  =  8.0,  where  the  G  value  is  defined 
as  the  number  of  particles  produced  per  100  eV  d-  posited. 

In  addition,  as  discussed  in  section  3.2,  each  positive  ion 
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created  will,  when  it  recombines  with  an  electron  or  negative 
ion,  produce  two  oxygen  atoms,  provided  that  other  reactions 
have  not  converted  the  ions  to  more  stable  species  like  0“ . 

O 

Since  G  (positive  ions)  =  3.25  (table  1),  ion  recombination 
can  increase  G(O^)  by  up  to  6.5,  giving  a  total  upper  limit 
of  G(03,  =  8.0  +  6.5  =  14.5. 

3.5  CHEMISTRY  COMPUTATION  METHOD. 

In  the  present  study,  time-dependent  chemistry  calculations 
have  been  made  for  situations  corresponding  to  published 
oxygen  radiolysis  experiments.  These  calculations  used  the 
CHEM  computer  code  developed  at  RDA,  which  employs  a  Gear 
integrator  (Ref.  51)  to  solve  numerically  the  coupled  differ¬ 
ential  equations  for  the  chemical  kinetics.  The  calculations 
employed  the  initial  production  rates  listed  in  table  1  and 
the  rate  coefficients  listed  in  table  2.  The  results  will  be 
discussed  in  the  next  section. 
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SECTION  4 


COMPUTED  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 

4.1  OZONE  PRODUCTION  IN  PULSE  RADIOLYSIS. 

Several  studies  have  been  carried  out  in  which  intense  short- 
duration  pulses  of  electrons  were  used  to  irradiate  oxygen 
(Refs.  1,  52-54).  The  electron  pulse  generators  used  had  peak 
electron  energies  of  600  keV  to  2  MeV  and  pulse  half-widths  of 
around  30  nsec . 

Typical  species  concentrations  calculated  for  the  pulse 
radiolysis  of  oxygen  at  1  atmosphere  are  shown  in  Figure  1 . 

The  principal  ions  are  0+ ,  0~  and  0~.  The  concentrations 
of  electrons  and  0 +  are  lower  due  to  efficient  attachment 
and  clustering  reactions.  The  ionic  charge  builds  up  during 
the  period  of  electron  irradiation,  reaches  a  maximum,  and 
then  decays  by  ion-ion  recombination,  approximately  as  t”1 
after  the  end  of  the  pulse.  At  late  times,  as  ozone  appears, 
the  ions  0^  and  0+  become  appreciable;  however,  by  this  time 
most  of  the  ionization  has  dissipated. 

Atomic  oxygen,  which  is  produced  directly  by  electron  bombard¬ 
ment,  accumulates  throughout  the  irradiation  period.  The 
deionization  processes  contribute  additional  atomic  oxygen, 
so  0  continues  to  build  up  early  in  the  deionization  phase. 

At  later  times  substantial  ozone  concentrations  are  produced 
from  the  association  of  0  with  •  Most  of  the  ionization 
has  disappeared  by  this  time.  The  molecular  oxygen  metastable 
states  and  decay  over  still  longer  time  scales,  and 
their  dissociating  reactions  with  0^  tend  to  maintain  a 
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CONCENTRATION  (molecules/cm 


TIME  (sec) 


Figure  1.  Time-dependent  species  concentrations  in  oxygen 
subject  to  a  brief  irradiation  pulse  varying 
with  time  as  t  exp(-t/t0)/  where  tQ  =  22  nsec. 
The  integrated  dose  is  1  x  10^  rads,  pressure 
760  torr,  and  temperature  293  K. 
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Figure  1.  Time-dependent  species  concentrations  in  oxygen 
subject  to  a  brief  irradiation  pulse  varying 
with  time  as  t  exp(-t/tQ),  where  tQ  =  22  nsec. 
The  integrated  dose  is  1  x  106  rads,  pressure 
760  torr,  and  temperature  293  K. 
concluded) . 
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low-level  concentration  of  atomic  oxygen,  but  this  has  little 
effect  on  the  0^  concentration  or  yield. 

In  table  3,  predicted  ozone  G  values  from  pulse  radiolysis  are 
compared  with  measured  values.  The  agreement  is  reasonably 
good,  considering  the  uncertainties  in  the  computations  and 
the  measurements.  The  principal  source  of  uncertainty  in  the 
calculations  is  probably  the  initial  yields  of  various  species 
in  oxygen  under  electron  bombardment.  All  of  the  pulse  radi¬ 
olysis  experiments,  and  the  corresponding  calculations,  show  a 
high  ozone  yield,  near  the  theoretical  limit.  Since  figure  1 
shows  that  at  atmospheric  pressure  most  of  the  charged-particle 
recombination  occurs  between  0+ ,  o~  and  0“,  this  high  ozone 
yield  implies  that  recombination  of  these  ions  yields  two 
oxygen  atoms,  as  assumed  in  Table  2.  Although  these  reactions 
are  not  listed  in  the  DNA  Reaction  Rate  Handbook  (Ref.  33) , 
the  similar  reaction,  0+  +  0“  -»  °2  +  °2 ,  is  included.  The 
present  results  strongly  suggest  that  the  products  should  be 
revised  to  0^  +  2  0. 

Similar  chemistry  calculations  were  made  for  pulse  radiolysis 

at  0.1  atm.  The  results  show  that  0+  ot  and  e  are  the  major 

z  4 

charged  species  in  this  case.  The  high  ozone  yield  observed 
at  this  pressure  by  Willis  et  al .  (Ref.  1)  then  indicates  that 
e  +  0+  yields  two  oxygen  atoms.  This  result  has  been  gener¬ 
ally  assumed  (Ref.  33)  but  not  previously  confirmed. 

4.2  EFFECT  OF  ELECTRON  SCAVENGERS  ON  OZONE  YIELD  IN  PULSE 
RADIOLYSIS . 


An  electron  scavenger  like  sulfur  hexafluoride  (SF^.)  can 
attach  electrons  to  form  SF“,  a  negative  ion  that  is  more 
stable  than  0“  or  0“,  and  hence  may  not  produce  oxygen  atoms 
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Table  3.  Ozone  yield  in  the  pulse  radiolysis  of  oxygen 
High  dose  rates:  1011  -  1013  rad/s  (1  rad  =  100  erg/g) 
Gas  temperature:  293  -  298  K. 


Type  of 

determination 

O2  pressure 

( torr ) 

G(03 ) 

raolec/ 100  eV ) 

Reference 

Measurement 

560 

13.0* 

52 

Measurement 

760 

13 . 8±0 . 7 

54 

Measurement 

30  -  800 

12 . 8+0 . 6 

1 

Calculation 

83,  830 

14.0+0.9 

1 

Calculation 

76,  760 

13.8-14.4 

This  work 

Theoretical  limit 

— 

14 . 5 

This  work 

* G  value  corrected  for  dosimetry  in  Ref.  1. 
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when  it  recombines  with  0+  or  0+.  These  and  other  relevant 

reactions  of  SFg  are  summarized  in  table  4.  Using  these  reac¬ 
tions  and  rate  constants,  chemistry  calculations  were  made 
for  the  pulse  radiolysis  of  oxygen  containing  1  percent 

(by  volume)  of  SF^f  corresponding  to  three  published  exper¬ 
iments  (Refs.  1,  52,  53) . 

Table  5  compares  the  measured  and  calculated  ozone  yields  . 

The  calculated  G(C>3)  =  8.0-8. 1  (varying  slightly  with  pres¬ 
sure)  lies  near  the  center  of  the  range  of  measured  values, 

G(0  )  =  6. 3-9.0.  Since  the  calculations  indicate  that  the 
major  recombining  ions  are  SF“ ,  0+  and  0  + ,  the  reasonable 
agreement  of  the  calculated  with  the  measured  ozone  yields 
confirms  the  assumption  (table  4)  that  recombination  of 
these  ions  does  not  yield  oxygen  atoms . 

4.3  OZONE  BUILDUP  DURING  CONTINUOUS  RADIOLYSIS. 

A  number  of  researchers  have  studied  the  effects  of  continuous 
radiolysis  of  pure  oxygen  (Refs.  1,  44,  55,  56).  Most  of 
these  studies  used  gamma  ray  sources  and  irradiation  intensi¬ 
ties  up  to  about  107  rad/hr,  with  total  doses  exceeding  10® 
rads  in  some  cases.  The  difference  between  the  energy  sources 
employed  in  the  pulse  and  steady  radiolysis  experiments  (elec¬ 
trons  versus  gamma  rays)  has  no  direct  bearing  on  the  ion 
chemistry  or  ozone  yield;  the  primary  ionization  products  per 
unit  energy  absorbed  are  essentially  the  same  in  both  cases. 
However,  the  steady  radiolysis  experiments  involved  much 
longer  irradiation  times  and  therefore  significant  ozone  accu¬ 
mulation  during  the  period  of  irradiation.  Consequently,  ions 
like  0~  and  0^  are  expected  to  be  more  important  here  than  in 
the  pulse  radiolysis  experiments . 
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Table  4.  SFg  chemistry. 


Reaction 

e  +  SF0  -♦  SFg 

SF6  +  02  **  SFg  +  02 

Rate  Coef f icient ( a) 

2.2  x  10-7 

7.0  x  10-11 

Reference  or 
Footnote 

57  ,  58 

34 

SFg  +  04**  SFg  +  2  02 

7.0  x  10-11 

b 

SFg  +  0  ■*  0“  +  SFg 

5.0  x  10-11 

34 

SFg  +  03  ->03  +  SFg 

3.2  x  10-11 

34 

SFg  +  0+  -*  SFg  +  0 

2.0  x  10-6 

b 

SFg  +  02  -*  SFg  +  02 

2.0  x  10-6 

b 

SFg  +  O4  ■*  SFg  +  2  02 

2.0  x  10-6 

b 

SFg  +  0£  -  SFg  +  02  +  03 

2.0  x  10-6 

b 

Footnotes : 


a.  All  rate  coefficients  have  units  of  cm^/sec . 

b.  Assumed  by  analogy. 
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Table  5.  Ozone  yield  in  pulse  radiolysis  of  oxygen  containing 
1  percent  (by  volume)  of  sulfur  hexafluoride;  temperatures  near 
295  K. 


Type  of  determination 

Pressure  (torr) 

G(03) 

References 

Measurement 

560 

* 

9 . 1 

52 

Measurement 

200-2000 

7 . 9 1 

53 

Measurement 

400,  700 

6 . 3+0 . 6 

1 

Calculation 

76,  760 

00 

o 

00 

h-> 

This  work 

Theoretical  limit, 
without  ionic  contrib. 

8.0 

This  work 

$ 

G  value  corrected  for  dosimetry  as  in  Ref.  1. 
tNormalized  using  G(O^)  =  12.8  in  pure  oxygen  at  760  torr 
( see  table  3 ) . 
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The  present  calculations  for  continuous  radiolysis  situations 
bear  out  this  expectation.  Although,  during  the  very  earliest 
period  of  the  irradiation  (corresponding  to  integrated  doses 

o 

of  less  than  10°  rads)  the  ozone  concentration  is  so  small 
that  little  0"  or  0  +  is  predicted  to  form,  by  the  time  the 
ozone  accumulates  to  levels  that  are  easily  measureable 
(around  104  rads)  0“  becomes  the  dominant  negative  ion,  and  0 
and  0+  the  dominant  positive  ions.  Consequently,  the  reason¬ 
able  agreement  of  the  calculated  and  measured  ozone  yields 
shown  in  table  6  supports  the  assumption  (table  2)  that  recom¬ 
bination  of  0*1  with  ol"  or  ot  does  not  lead  to  dissociation, 
o  4  5 

(The  discrepancy  of  about  25  percent  between  the  calculations 
and  the  measurements  is  not  readily  explainable  but  probably 
has  little  practical  significance.) 

As  the  oxygen  irradiation  time  is  lengthened  and  the  total 
dose  increases  above  about  10  rads,  the  calculated  ozone 
production  rate  slows,  due  primarily  to  the  assumed  process, 

O3  +  O3  -*  O2  +  ^  (Ref.  43),  which  effectively  destroys 
two  ozone  molecules  per  reaction.  (The  well-known  ozone 

destruction  reaction,  0  +  0^  -♦  2  O2 ,  is  much  slower  under 
most  radiolysis  conditions.)  A  qualitatively  similar  behavior 
has  been  observed  in  two  sets  of  high-dose  radiolysis  measure¬ 
ments  (Ref.  44,  56).  However,  Kircher  et  al .  (Ref.  44) 

used  "commercial"  oxygen  containing  impurities  that  may  have 
affected  the  ozone  chemistry;  moreover,  their  pressures  were 
much  higher  than  atmospheric.  Consequently,  no  comparisons 
of  Kircher' s  measurements  with  chemical  calculations  have 
been  attempted. 

Sears  and  Sutherland  (Ref.  56)  carried  out  experiments  near 
atmospheric  pressure  using  both  all-glass  vessels  and  quartz 
vessels  having  stopcocks  lubricated  with  fluorocarbon  grease. 
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Table  6.  "Initial"  ozone  yield  in  the  steady  r^diolysis  of 
oxygen.  ("Initial"  implies  low  total  dose,  ~104  rads.) 

Gas  temperatures  are  near  295  K  unless  otherwise  indicated. 


Type  of 
determination 

02  pressure 

g(°3) 

References 

Measurement 

( 2- 10 ) xlO4 

* 

5 . 9 

44 

Measurement 

(200  K) 

760 

* 

6.7  t 

55 

Measurement 

665 

6t* 

56 

Measurement 

(77  K) 

not  stated 

6. 2+0 . 6 

1 

Theoretical 

estimate 

— 

7 . 5+0 . 9 

1 

Calculation 

665 

8-9 

This  work 

Theoretical  limit, 
without  ionic  contrib. 

8 . 0 

This  work 

G  value  corrected  for  ozone-iodide  stoichiometry  in  Ref.  1. 
tG  value  corrected  for  dosimetry  in  Ref.  1. 

*See  Ref .  1 . 
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They  observed  that  the  ozone  yields  were  greater  in  the  ves¬ 
sels  with  stopcocks.  This  observation  is  consistent  with  the 
fact  that  many  fluorine  compounds  will  remove  the  negative 
charge  from  ions  like  C>3 ,  and  hence  interfere  with  the  C>3  +  03 
destruction  mechanism. 

Unfortunately,  Sears  and  Sutherland  reported  dose-dependent 
measurements  only  for  the  vessels  with  greased  stopcocks. 

Their  results  are  plotted  in  figure  2,  along  with  calcula¬ 
tions  for  pure  oxygen  and  oxygen  contaminated  with  0.0044 
percent  (44  ppmv)  of  a  hypothetical  electron  scavenger. 

The  scavenger  is  assumed  to  have  the  same  reactions  as  SFg 
(table  4),  except  that  the  charge  transfer  reactions  with  0~ 
and  03  proceed  in  the  opposite  direction  (i.e.,  the  scavenger 
has  a  greater  electron  affinity  than  0  and  03) .  The  experi¬ 
mental  points  in  figure  2  fall  between  the  calculated  points 
for  pure  and  contaminated  oxygen,  strongly  suggesting  that  the 
stopcock  grease  is  producing  an  electron  scavenger  that 
removes  the  charge  from  0- 

4.4  STEADY- STATE  OZONE  CONCENTRATIONS  IN  CONTINUOUS 
RADIOLYSIS . 


With  long  enough  irradiation  times  (corresponding  to  doses 

7 

above  approximately  10  rads)  the  chemical  processes  reach 
a  steady  state  where  ozone  is  destroyed  as  fast  as  it  is 
formed.  The  ozone  concentration  thus  eventually  approaches 
a  constant  value.  Sears  and  Sutherland  (Ref.  56)  present 
measurements  for  this  concentration  in  all-glass  (uncontam¬ 
inated)  vessels,  as  a  function  of  the  dose  rate.  These 
values  are  compared  in  figure  3  with  values  calculated  using 
the  reaction  rates  in  table  2.  The  agreement  is  fairly 
good.  The  experimental  values  vary  approximately--and  the 
calculated  values  almost  exactly--wi th  the  square  root  of 
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OZONE  CONCENTRATION  (ppmv) 


DOSE  (rads) 

Figure  2.  Measured  and  calculated  ozone  concentrations 
versus  cumulative  radiation  dose  for  the  con¬ 
tinuous  (slow)  radiolysis  of  oxygen  at  665  torr 
and  283  K.  The  diagonal  lines  correspond  to  loci 
of  constant  G  values.  For  the  calculated  points 
connected  by  vertical  dashed  lines,  the  lower 
values  correspond  to  pure  oxygen  while  the  upper 
values  correspond  to  oxygen  contaminated  with 
44  ppmv  of  a  hypothetical  electron  scavenger 
(see  text) . 
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STEADY  STATE  OZONE  CONCENTRATION  (ppmv) 


DOSE  RATE,  6  (rads/hr) 


Figure  3.  Comparison  of  measured  and  calculated  steady-state 
ozone  concentrations  in  irradiated  oxygen  as  a 
function  of  dose  rate. 
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the  dose  rate  (solid  line  on  figure  3).  As  pointed  out  by 
Sears  and  Sutherland,  this  behavior  expected  when  the  pri¬ 
mary  removal  process  is  the  reaction  of  a  major  ion  with  0 

O  f 

since  in  a  steady  state  dominated  by  ion-ion  recombination 
the  ion  concentration  (of  either  sign)  is  proportional  to 
the  square  root  of  the  dose  rate . 

The  agreement  of  the  steady-state  measurements  with  the 
calculations  tends  to  support  the  assumed  rate  coefficient 
and  products  of  the  reaction  0"  +  o,  oi  +  2  0„  (table  2) 
However,  before  this  support  can  be  considered  strong,  alter¬ 
native  destruction  reactions  must  be  ruled  out.  Note  first, 
that  to  fit  the  radiolysis  observations,  the  species  reacting 
with  and  destroying  ozone  must  be  an  odd  oxygen  species,  like 
Og,  0+  or  0(*D);  excited  molecules  can  only  break  up  ozone 
to  give  0  +  C^,  where  the  0  atoms  immediately  re-form  ozone. 
The  decrease  in  ozone  destruction  in  the  presence  of  fluoro¬ 
carbon  grease  (figure  2)  strongly  suggests  a  negative  ion, 
since  many  fluorine  compounds  have  large  electron  affinities 
and  can  scavenge  charge  from  other  negative  ions.  Moreover, 
if  a  neutral  species  (other  than  ozone)  were  the  reactant,  its 
primary  removal  process  would  have  to  be  reaction  with  itself 
or  with  another  free  radical  or  ion,  in  order  to  yield  the 
observed  square  root  behavior.  Enough  is  known  about  the 
potentially  significant  neutral  species  (in  fact,  they  have 
been  included  in  table  2  and  in  the  calculations)  to  show  that 
none  satisfy  this  requirement . 

Possible  ozone  destruction  by  ions  like  0-,  0+  and  0^  should 
also  be  considered.  However,  the  calculations  show  that  the 
concentration  of  0“  in  the  radiolysis  experiments  is  always 
much  too  small  to  contribute  significantly  to  the  ozone  loss. 
The  0+  concentration  is  also  small,  because  it  is  formed  only 


38 


from  0  +  ,  which  has  a  small  concentration,  and  it  is  rapidly 
destroyed  by  charge  transfer  to  02-  (For  these  reasons,  O3 
was  not  included  in  the  present  chemical  calculations.) 

On  the  other  hand,  the  concentration  of  0  +  is  large  enough 
for  it  to  destroy  significant  amounts  of  ozone,  if  the  rate 
coefficient  were  large  enough.  At  the  lowest  and  highest 
dose  rates  employed  by  Sears  and  Sutherland,  0+  makes  up  about 

O 

20  and  60  percent,  respectively,  of  the  total  positive  ions, 
the  remainder  being  0*,  which  cannot  destroy  ozone.  Hence,  if 

the  0+  +  03  and  O3  +  03  coefficients  were  equal,  and  both  were 
a  factor  1/1.6  smaller  than  the  value  used  for  the  latter  reac¬ 
tion  in  the  calculations  shown  in  figure  3,  the  low  end  of  the 
curve  would  be  raised  about  a  factor  of  1.6/1. 2,  while  the  high 
end  would  remain  unchanged  (neglecting  the  fairly  small  effect 
of  the  changed  ozone  amount  on  the  0+/0+  ratio).  This  modified 
prediction  would  still  fit  the  measurements  reasonably  well. 

However,  if  the  0+  rate  coefficient  were  twice  the  0“  coeffi- 

o  3 

cient,  there  would  be  a  significant  discrepancy. 

It  follows  that  the  0  +  +  0^  rate  coefficient  is  probably  less 
than  6  x  10  cm  /s  and  may  be  much  smaller.  In  addition,  the 
value  used  for  the  0~  reaction  in  the  calculations  1  x  lO-*2 

q 

cmJ/s,  is  probably  an  upper  limit,  since  it  would  have  to  be 
lowered  if  0+  contributes  significantly  to  the  destruction, 
or  if  the  ion-ion  recombination  coefficient  is  smaller  than  the 
value  of  2  x  10-h  cm3/s  used  in  the  calculations  (compare  Ref.  50) , 
since  this  would  increase  the  0“  concentration.  Accordingly, 
a  reasonable  choice  for  the  room-temperature  rate  coefficient 
of  the  0 '  +  c>3  reaction  is  5  x  10  cm°/s ,  which  may  be  correct 
to  within  a  factor  of  2 .  It  probably  varies  with  temperature 
roughly  as  exp( -1000/T ) ,  as  deduced  by  Kircher  et  al.  (Ref.  44) 
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for  the  temperature  dependence  of  ozone  destruction  in  high- 
pressure  oxygen  radiolysis.  However,  since  the  temperature  may 
also  shift  the  relative  ion  concentrations  and  change  the  ion- 
ion  recombination  coefficient,  the  observed  temperature  depen¬ 
dence  may  not  be  attributable  solely  to  the  0~  +  o  reaction. 

o  3 
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SECTION  5 


CONCLUSIONS 

The  detailed  electron-impact  and  chemical  calculations 
carried  out  in  this  study  have  given  results  which  agree 
within  25  percent  with  published  measurements  of  ozone 
production  in  a  variety  of  radiolysis  experiments.  The 
calculations  employed  electron-impact  cross  sections  newly 
deduced  from  published  data,  and  standard  sets  of  reaction 
rate  coef f icients ,  with  the  addition  of  the  following  ion- 


ion  recombination  reactions: 

Oj  +  0"  -  2  0  +  2  02,  (18) 
°4  +  °4  -  2  0  +  3  02,  (19) 
°4  +  °3  ■*  2  °2  +  °3  '  (20) 
°5  +  °3  -  °2  +  2  °3-  (21) 


(Some  other  recombination  reactions  were  also  added  for 
completeness,  but  they  have  a  negligible  effect  on  the 
results . ) 

In  conventional  ion  recombination  experiments  it  has  not 
been  possible  to  identify  the  products  of  these  reactions. 
The  present  analysis  of  the  radiolysis  experiments  permits 
confident  assignment  of  these  products  (except  that  in  the 
last  two  reactions  the  might  be  replaced  by  0  +  0o) . 
Moreover,  the  products  of  the  first  reaction  strongly  imply 
that  the  0+  +  0”  recombination  reaction  similarly  yields 
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2  0  +  °2»  instead  of  2  02  as  conventionally  assumed.  In 
addition,  the  analysis  of  the  steady-state  ozone  concen¬ 
tration  attained  by  prolonged  irradiation,  as  a  function 
of  the  dose  rate,  suggests  that  the  ozone-destroying  reaction 

+  Og  -»  o2  +  2  (15) 

occurs  with  a  rate  coefficient  within  a  factor  of  two  of 
5  x  10-13  cm3/s  at  room  temperature,  and  probably  with  a 
temperature  variation  of  about  exp  (-1000/T).  The  reaction, 

°5  +  °3  ■*  °4  +  2  °2  (22) 

may  also  provide  a  significant  ozone  sink.  However,  it 
cannot  be  the  dominant  loss  mechanism,  so  its  rate  constant 
cannot  be  larger  than  about  6  x  10  10  cm  /s. 

Calculations  for  oxygen  containing  traces  of  electron  sca¬ 
vengers  (such  as  fluorine  compounds)  show,  in  agreement  with 
experiment,  that  these  scavengers  can  affect  the  ozone  yield 
in  two  opposing  ways:  In  pulse  radiolysis  they  can  reduce 
the  ozone  yield  by  changing  the  ions  to  species  that  produce 
less  atomic  oxygen  when  they  recombine.  In  continuous  irra¬ 
diation,  where  the  ozone  builds  up  and  itself  reduces  the 
ozone  yield  by  electron  scavenging,  addition  of  a  stronger 
scavenger  can  remove  the  charge  from  0”  thus  inhibiting 

O 

the  O3  +  O3  reaction  that  destroys  ozone,  and  increasing 
the  ozone  yield. 
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